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Abstract. A method for gaining information about the phonon-number 
moments and the generalized nonlinear and linear quadratures in the motion 
of trapped ions (in particular, position and momentum) is proposed, valid inside 
and outside the Lamb-Dicke regime. It is based on the measurement of first time 
derivatives of electronic populations, evaluated at the motion-probe interaction 
time T = 0. In contrast to other state-reconstruction proposals, based on 
measuring Rabi oscillations or dispersive interactions, the present scheme can 
be performed resonantly at infinitesimal short motion-probe interaction times, 
remaining thus insensitive to decoherence processes. 
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1. Introduction 

A single ion in a radio-frequency trap is an ideal system for the investigation of 
the basic and intriguing features of quantum mechanics J^. Recent advances in 
the manipulation of the internal and external degrees of freedom of the trapped 
particle allowed, for example, the realization of laser cooling to the ground state 
for extended periods of time [IIISIt the creation of nonclassical motional states 01 El, 
entanglement between external and internal degrees of freedom of one ion or several 
ions 13 El m ' quantum gates for quantum computation El E] , and single 
photon emission on demand |14| . 

An important aspect at this advanced level of quantum state engineering is the 
ability to obtain as much information as possible, and with high efficiency, about the 
motional state of the trapped atom. Moreover, since the Lamb-Dicke (LD) regime, 
where the motion of the particle is restricted to a region smaller than a wavelength, is 
sometimes difficult to achieve in a trap, it is desirable to have measurement schemes 
working also outside this regime |15|. 
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So far, several techniques have been proposed to gain information about the 
motional quantum state of a trapped particle. The central idea in all these methods 
is the mapping of the ionic external dynamics to an internal degree of freedom, where 
the latter can be read out, for example, with electron shelving techniques ^Hl- Some 
of these procedures are applicable only inside the LD regime ^1 El El HOI i whereas 
others work also beyond this limit [HI E21 ESI ElESl ES| , allowing to derive m some 
cases the complete motional density operator. However, a typical requirement of 
conventional schemes is that the measurements have to be performed during relatively 
long interaction times, through several Rabi oscillations in the resonant case or slow 
phase shifts in the dispersive case, where decoherence mechanisms cannot be neglected. 

In this paper, we present a method that allows to determine the expectation value 
of the phonon-number moments (n^) of a trapped ion for any positive integer p, as well 
as any generalized motional nonlinear quadrature ^{f{h)ae~^'^ + f{n)e'^'^) \. Here, 
a and a) are the phonon annihilation and creation operators, respectively, n = a^a, 
/(n) is a function of n, and (j) is an arbitrary phase factor. In particular, we are 
able to measure the expectation value of the linear quadratures, i.e., when /(n) — 1, 
including the ionic position and momentum. Our technique is not restricted to the LD 
regime and, moreover, can be extended easily to all spatial dimensions, or to systems 
containing more than one trapped ion. Furthermore, in contrast to the methods 
presented so far, our proposal allows to obtain the required information in an extremely 
short motion-probe interaction time, which is useful in particular in presence of strong 
decoherence processes. 

In Sec. 121 we show for a single trapped ion how to obtain information about the 
expectation values of the phonon-number moments (Sec. I2.1|l and the generalized 
motional nonlinear and linear quadratures (Sec. I2.2|l through infinitesimal-time 
motion-probe interactions. In Sec. we extend the method to the TV-ion case. In 
Sec. 01 we summarize our central results. 

2. Single ion case 

2.1. Phonon-number moments 

We consider a single two-level ion trapped in a harmonic potential. For the sake of 
simplicity, we confine our treatment to one motional degree of freedom, though all 
results are easily generalizable to three dimensions. The system is described by the 
Hamiltonian 

HQ = hv{h + ll'2.) + hujQ\e){e\, (1) 

where v is the harmonic oscillator frequency, |e) and \g) are the electronic upper 
and lower states of the two-level ion, respectively, and o^o is the associated transition 
frequency. 

The ion is excited by a travelling laser beam resonant with the carrier electronic 
transition, leaving unchanged the motional populations. In the usual rotating-wave 
approximation, and in the interaction picture, the Hamiltonian reads |15l I24| 

i?int = \hnL{^+ + <T-)/o(n; 77), (2) 

I Generalized nonlinear quadratures can be related to nonlinear coherent states, theoretically studied 
in 1271 . A method for generating nonlinear coherent states in trapped ions, but not for measuring 
them, can be found in 1281 . 
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where fi^ is the Rabi frequency, (t* are the electronic two- level flip operators, ry is 
the LD parameter (77 = fcj.(a;o), (xg) being the extension of the ground state of the 
motional mode and the projection of the laser wave vector on the trap axis), and 
/o(n;7?) is given by 



fo{n;v) = e 



[n~l)V 



(3) 



by 



At any time t, the probability of finding the ion in the excited state |e) is given 



PeW=Tr p(t)\e){e\ 



(4) 



where p{t) is the system density operator, describing the internal and external degrees 
of freedom of the trapped particle. 

Considering that, for any operator A, 



±IA)^U\A.H 
at in 



\ dt r 



\e){elH 



we get, for A = |e)(e|, 

dPe ^ 1 

dt ih 

In our case, we have [|e)(e|,i?o] — 0, so that 

dPe _ 1 
dt ih 
Since, in the interaction picture 



[|e)(e|,^int] 
we obtain, from iQ and ©, 



e){e\,Hint 

((7+ - (7")/o(ri,?7), 



dPe ^ ^L, 

dt ~ 2i 



= ^Tr[p(i)(a+-a-)/o(n,?7)]. 
Next, we consider that the ion is prepared initially in the state 

m = \±4,){±^\(g,pf, 

where 

\±^) = ^{\g)±e^*\e)), 

and Pf is the phonon state we aim to characterize. Then, a straightforward calculation 
yields 

Tr[p(0)(a+-a-)/o(n,ry)] = T^sin(</>)Tr[p//o(n, r/)] , (12) 
and thus, using Eq. 

dP^ 



(5) 

(6) 

(7) 

(8) 

(9) 
(10) 

(11) 



dr 



= Tsin(0)(/o(?T.,?7)), 



(13) 



T = 



where t is the dimensionless time VLLt/2 and the token ±0 stands for the two 
parameters (sign ± and phase 0) defining state ji^) of Eq. Hll|) . 
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Eq. shows that the mean value of the nonhnear operator fo{n',ri) is 

determined by the time derivative of the population of the excited state at the initial 
interaction time r = 0. This allows to gain information about the ionic motional state, 
in particular, as will be shown below, about the phonon-number moments. Moreover, 
since the time derivative of the excitation probability Pe * is evaluated at r = 0, this 
information can be obtained in a very short motion-probe interaction time, even before 
decoherence mechanisms can affect the initial coherent evolution. Remark further that 
the "contrast" in the measurement of {fo{n, rj)), as seen in Eq. can be tuned with 
the proper choice of the phase (f>. 

For small LD parameters, a series expansion of the nonlinear operator fo(ji',ri) 
leads to the expression 

{foih;v))^l-v'{n) + ^{n'). (14) 

Thus, if we repeat the measurements with two known LD parameters, rji and 772 
(varying, e.g., the angle between the laser beam and the trap axis), we can derive (h) 
and (fp) by solving the linear system 

(15) 

from which, for example, the Fano-Mandel parameter |29| 

(n) 

can be extracted. It is known that the Fano-Mandel parameter allows to determine 
the "classicality" of a given motional state, since quantum states with Q > 1 are 
considered as classical states, and those with Q < 1 as nonclassical states, due to their 
sub-Poissonian phonon-number distribution. 

This low order moment determination scheme can be generalized to higher order 
vibronic moments (nP) {p > 0). Indeed, for higher LD parameters, the series expansion 
of fo{n, rj) in Eq. Q will extend beyond the {tt?) term, up to a certain moment {fi^) §. 
We could repeat then the measurement technique, as in Eq. (|15() . for N known LD 
parameters, yielding a linear system from which the N first vibronic moments become 
accessible. Each (n^), < p < N, could also be determined by measuring only once 
the initial time derivative of P^*' . To show this, we require that the ion is submitted to 
N simultaneous laser interactions, each of them resonant with the electronic transition 
and leaving the motional state unchanged. In the interaction picture, they give rise 
to the simultaneous action of Hamiltonians 

H',^,^^nn,a+fo{n;Vj)+^-c., j^l,...,N, (17) 

where flj and rjj are, respectively, the electronic Rabi frequency and the LD parameter 
of laser j. The total Hamiltonian will be given in this case by 

Hint = E Ant = ^(a+ +a-)Fo(n), (18) 
J = l 

§ The series of Eq. convergent for all rj. When the sum of this series is truncated at Z = A'^, 

the coefficient of the term is given by {—l)^e~^ ^'^rj^^ /N^^ . For any iq value, this coefficient 
can be made arbitrarily small by a proper choice of and the series may be truncated within good 
approximation up to this term. 
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where fl^ = max(rij) and 



Fo{n) 



N 



/^7^/o(";»7j) 



If the system is again initially prepared in the state ji^) 
to Eq. (O, 



(19) 

D/5/, we get, similarly 



dPe 



dr 



Tsin(0)(Fo(r^)). 



It has been shown by de Matos Filho and Vogel 
in the form of a Taylor series 



Fo{n) 



OO 

p=0 



(20) 

that F(){n) may be rewritten 
(21) 



where the Taylor coefficients Cp are given by 



c„ = < 



N 



(-i)^E' 



^ o 

J = l 



-Tjf /2 % 



„2m 



il2 



if 



if p^O, 



(22) 



with 



m— 1 

E 

,<...<x 



JilJi 



■3i 



if 



if p < m. 



(23) 



This means that for given values of the N LD parameters 77^ (fixed by the geometry 
of the laser beams regarding to the trap axis up to the maximum value ?7max = k{xo), 
k being the wave vector of the lasers) , the Cp coefficients are linear combination of all 
Rabi frequencies flj . In this way, the use of N lasers allows to fix at will N coefficients 
of the Taylor series. In particular, if the N first coefficients are set to 0, we obtain 



Fo{h) = Oi 



max 



)• 



(24) 



2N 

Note that it is always possible to choose the N value such that 0{^^^fffrn^) is negligible, 
even outside the LD regime where ?7max can be greater than 1. The mean value of 



-AT 



) can be verified experimentally by measuring the time derivative of the 
population of the excited state at interaction time r = 0. According to Eq. H2U|) . 

2N 

the outcome of this measurement yields (©(^^p-n^)) so that we can check if it is 
indeed negligible for the phonon distribution we aim to characterize (in this case it is 
recommendable to set (p = ±7r/2). 

The Rabi frequencies of the N lasers could also be chosen in such a manner that 
only a single coefficient Cp {p < N) is equal to 1, the others remaining equal to 0. 
In this case, according to Eq. I|21(l and considering that all terms h^, r ^ N, are 
negligible. 



Foin) 



(25) 
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± 

and the measurement of 



yields directly the vibronic moment {fiP). This 

T = Q 

step can be reproduced for any p < N. 

It is noteworthy to mention that the knowledge of the phonon-number moments 
(tV) for all positive integers p allows to derive the complete phonon distribution 
p{n) = {n\pf\n), as known from classical statistics Also, the generalization to the 
31? case provides measurement schemes for (n^), (%), {nz}, {n"^}, (^y), ("-z)j {nxfi-y), 
{fixfiz), (hyhz), (hxhyhz), and so on. Recently, the relevance of measuring different 
photon-number moments, through quantum field homodyning, for determining the 
"classicality" of arbitrary quantum states has been considered . 

2.2. Generalized nonlinear quadratures 

Next, we show that by using a red or blue sideband excitation, expectation values 
of generalized nonhnear quadrature moments, ^{f{h)ae~^'^ + a^/(n)e*'^), can be 
determined using similar techniques. In particular, we have access to the expectation 
values of the generalized linear quadrature moments, among them position (x) and 
momentum (p) \\. As will be shown below, these expectation values can be obtained 
inside and outside of the LD regime. 

Let us consider that the ion is excited with a red-sideband detuned laser (the blue- 
sideband case can be treated equally, leading to similar results). In the usual rotating- 
wave approximation, and in the interaction picture, the Hamiltonian reads |15[ I24| 

i?mt = ^fif^Lf?<T+/i(n;77)a-f H.c, (26) 

where 

If the ion system is initially prepared in the state |±0)(±0| ® pf, we obtain 
straightforwardly from Eq. © 



dPt' 



dr 



±i(A(n;r;)ae-'* + at/i(";r/)e**), (28) 



where r is the dimensionless time riflLt/2. In the LD regime, /i(n; 77) = 1 and Eq 
reduces to the generahzed quadrature = (ae"*"^ -I- a)e'^'^)/2, 

dPt 



dr 



__^ = MX,) = ±l{cosW^^+M^)^), (29) 

where x and p are the trapped ion position and momentum, respectively, and (xq) 
and {po) are the spreads of these quantities in the ground state of the trap potential. 
Eq. (|29|l shows that any quadrature moment can be experimentally determined 
from the measurement of the initial time derivative of the population of the excited 
state, whereas, in particular, (x) and (p) can be obtained by choosing </> = 0, tt and 
(f) = ±7r/2, respectively. 

II In this limit our results for trapped ions are analogous to the ones derived for measuring field 
quadratures in cavity QED setups: see 1321 . 
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Outside the LD regime, a similar information is gained using N simultaneous red- 
sideband detuned lasers. In the interaction picture, they give rise to the Hamiltonians 

Ant = ',^m,r^jij+ h{n; + H.c, j = 1, . . . , TV, (30) 

where Q,j and r]j are, respectively, the electronic Rabi frequency and the LD parameter 
of laser j. The total Hamiltonian is given in this case by 

N 

Hint = J2 Ant = lhnLa+Fi{h)a + H.c, (31) 



2 



where flL = max(rij) and 



N 



If the ion is initially in the state |±0)(±0| (8) p/, we obtain, similarly to Eq. H28|) . 
with T = nLt/2, 

±i (^Fi{n)ae-"'' + a) Fi{n)e"l''^ . (33) 



dPe 



dr 

According to de Matos Filho and Vogel 24 , i^i (n) may be again rewritten in the 
form of a Taylor series 

oc 

A(ri) = ^CpnP, (34) 

p=0 

where the Cp coefficients are linear combination of the N laser Rabi frequencies flj . 
Similarly to the case of a resonant laser described in the preceding section, we can 
choose the TV value in such a way that 0{cNn^) is negligible and fix next the TV 
Rabi frequencies to set all coefficients Cp (p < TV) to 0, except co to 1. In this case, 
Fi{n) ~ 1, Eq. (|33|l reduces to Eq. (|29|l and the quadrature moments (and more 
specifically (x) and (p)) are determined as done in the case of the LD regime. Clearly, 
by following a similar procedure, we could also engineer Fi{h) to describe an arbitrary 
polynomial. 



3. TV-ion case 



We now briefly discuss how our proposal can be generalized to a chain of TV identical 
two-level ions in a linear Paul trap. This system is described by the Hamiltonian 

N N 

Ho = Y, + 1/2) + fi^o 51 Ie)fc(e|> (35) 

j=i k=i 

where i^j are the frequencies associated with the collective motional modes, iij = a^dj 
are the respective phonon-number operators, \e)k and \g)k are the electronic states of 
the two- level ion k, and loq is the electronic transition frequency of each ion. 

The ions are illuminated by a laser beam resonant with their electronic 
transition while leaving unchanged the phonon population, realizing a nonlinear carrier 
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excitation. In the usual rotating-wave approximation and in the interaction picture, 
the Hamihonian reads ^1 1^ 

TV N 

Hint = 2 E ('^fc /o {n, ;Vj), (36) 

k=l j=l 

where VIl is the electronic Rabi frequency, are the electronic two- level flip operators 
of atom fc, r]j is the LD parameter related to the collective mode j {rjj — k,j;{xo)j, (xo)j 
being the extension of the ground state of mode j and the projection of the laser 
wave vector on the trap axis), and 

In the following, we denote the product IljLi /o(%S ^j) by J^q. The generalized 
relation 



diPe)k 1 



\e)k{elH 



(38) 



dt ih 

can be easily obtained, as in Eq. (|SJ), where (Pe)fc is the probability of finding ion k 
in its excited state. As |e)fe(e| commutes with Hq and any operator associated with 
other ions fc', we get immediately 

diPe)k ^ l_ 

dt ih 

and thus 

d{Pe)k _ 



(^[\e)k{el^hnUai +a^)To]), (39) 



dt 2^^'^P^*^^'^^''''^^^''^' ^^^^ 

where p{t) is the iV-ion system density operator. 
Let us consider the following initial state 

p{0) ^ Pk<» PA<» Pf, (41) 

where 

Pk = \±4.)k{±4>\ (42) 
is the electronic density operator of the fc-th ion, with 

|±^)fc = i=(|5)fc±e^1e),), (43) 

Pa is an arbitrary electronic density operator of the remaining ions, and Pf is the 
collective motional density operator associated with the N eigenmodes. Following 
similar steps as in the previous sections, we obtain 

Tr[p(0)(a+ ~ a^):Fo] = sin(0)Tr [p/./-o] , (44) 

similar to Eq. ifT^ . and thus, using Eq. (^nj, 

d{Pc)k 



dr 



= Tsin(</.)(.Fo), (45) 

T = 

where r is the dimensionless time t = Q.i,t/2. Eq. H45|) tells us that by measuring the 
time derivative of the population of the excited state of one ion at t = 0, we are able 
to gain information about the collective nonlinear operator Tq. 
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4. Summary 

In conclusion, we have proposed a method that ahows to determine the phonon- 
number moments and the motional nonlinear and linear quadratures of a trapped 
ion, in particular the ionic position and momentum. Our method makes use of the 
nonlinear behavior of the ion-laser interaction in harmonic traps. The measurement 
of the phonon-number moments (n^) requires resonant carrier interaction, with no 
phonon gain or loss, while the measurement of the nonlinear quadratiirc moments 
i(/(n)ae-^'^ + a)f{n)e^'t') demands the use of red or blue sideband excitations. In 
contrast to methods presented so far, our proposal is designed for measurements 
realized in very short probe laser interaction times, thus preventing the noisy action 
of decoherence processes. In addition, we have shown that our scheme works inside 
and outside the LD regime, and that it can be generalized to the case of N ions. In 
the latter case, information about a collective property is gained by the measurement 
of the excited state population of a single ion. 
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